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Abstract 

To find the Hermitian phase operator of a single-mode electromagnetic field in quan- 
tum mechanics, the Schrodinger representation is extended to a larger Hilbert space 
augmented by states with infinite excitation by nonstandard analysis. The Hermitian 
phase operator is shown to exist on the extended Hilbert space. This operator is natu- 
rally considered as the controversial limit of the approximate phase operators on finite 
dimensional spaces proposed by Pegg and Barnett. The spectral measure of this oper- 
ator is a Naimark extension of the optimal probability operator- valued measure for the 
phase parameter found by Helstrom. Eventually, the two promising approaches to the 
statistics of the phase in quantum mechanics is synthesized by means of the Hermitian 
phase operator in the macroscopic extension of the Schrodinger representation. 

1. Introduction 

The existence and properties of a Hermitian operator on a Hilbert space corresponding to 
the phase of the electromagnetic field has provoked many discussions since Dirac first 
discussed the problem. According to the uniqueness theorem of the irreducible representa- 
tions of the canonical commutation relation due to von Neumann, the commutation relation 
between the number operator and the phase operator which Dirac presupposed from the cor- 
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respondence between the commutator and the classical Poisson brackets cannot be realized. 
Further, Susskind and Glogower |2| clearly demonstrated that the polar decomposition of 
the annihilation operator into the unitary operator of the exponential of the phase and the 
square root of the number operator presupposed by Dirac is also impossible. However, Pegg 
and Barnett recently made an interesting proposal for the problem. They constructed 

approximate Hermitian phase operators on finite dimensional spaces and claimed that the 
statistics obtained by their operator approaches the statistics of the phase as the dimension 
tends to infinity. However, they have failed to find the Hermitian phase operator on an 
infinite dimensional space as the limit of their approximate operators. 

On the other hand, another approach to the problem has been established in quantum 
estimation theory || This theory discusses optimization problems of quantum measure- 
ments quite generally. The statistics of measurement is represented in this theory by a 
probability operator-valued measure (POM) on a Hilbert space which extends the conven- 
tional description by a Hermitian operator. In this approach, the optimum POM of the 
estimation problem of the phase parameter was found by Helstrom ||, and mathematically 
rigorous development of this approach is given by Holevo [[7]. 

A promising aspect shared by these two approaches is that the statistics of the phase 
obtained by the limit process of Pegg and Barnett coincides with the one represented by the 



optimum POM of the phase parameter PI— 1|12|| . This shows, however, that contrary to their 



claim the limit of the exponentials of the approximate phase operators is nothing but the 
well-known Susskind-Glogower exponential phase operator ||, as long as the limit is taken 
on the Hilbert space of quantum states with the weak operator topology. The limit in the 
weak operator topology does not preserve the product operation and hence demolishes the 
desired properties of the limit operator which Pegg and Barnett || described intuitively. 
According to the Naimark theorem, every POM can be extended to a projection- valued 
measure on a larger Hilbert space which gives rise to a Hermitian operator by the spectral 



theory representing an observable in the standard formulation of quantum mechanics. This 
suggests that there exists the Hermitian phase operator somewhere beyond the Hilbert space 
of quantum states. Thus in order to realize the intuitive limit of the approximate phase 
operators, we need an alternative mathematical construction other than the limit on a Hilbert 
space. 



In this paper the attempt from nonstandard analysis outlined in |L2[] is developed for 



this purpose. The nonstandard analysis was invented by Robinson JT3[ and has yielded 
rigorous and fruitful mathematics of infinite and infinitesimal numbers. We construct a 
natural extension of the Schrodinger representation and show that the desired Hermitian 
phase operator exists on this extended Hilbert space. The Hilbert space of this extension of 
the Schrodinger representation is the direct sum of the original space of quantum states and 
the space of states with infinite excitation which are naturally considered as the classical 
limits of the ordinary quantum states. 

In the conventional approach, microscopic properties and macroscopic properties are 
discussed separately in quantum mechanics and in classical mechanics. Although the corre- 
spondence principle bridges both mechanics by the mathematical process of taking the limit, 
this approach cannot describe the quantum mechanical coherence between microscopic states 
and macroscopic states. Our new representation realizes such a coherent description of quan- 
tum and classical mechanics; it is in such a representation that the phase operator behaves 
as a Hermitian operator. Obviously, the present method is applicable to other difficulties in 
quantum mechanics concerning the bounded dimensionless quantities such as the rotation 
angle and the quantities associated with them such as the time of periodic motions. A po- 
tential application of this representation other than the above problems is the measurement 
problem, where the unitary time evolution in an amplifier evolves from a quantum state to 
a state with infinite excitation Applications to these problems will be discussed in the 



forthcoming papers. 



For bibliography on the phase operator problem we shall refer to the references of p|, O, 



KJ, and for the recent developments emerged by the Pegg-Barnett proposal the references of 



[I7|. For quantum estimation theory [El |7|, for quantum measurement theory [JTS , B5|, and 



for operator algebras PU[ . For basic methods of nonstandard analysis, we shall refer to ETIJ . 



Applications of nonstandard analysis to physics is not new and has been developed in such 
papers as and in monograph 



2. The Susskind-Glogower operators 

The single-mode electromagnetic field is a well-known physical system which is modeled by 
the quantum mechanical harmonic oscillator with unit mass. Let Ti, be the Hilbert space of 
the Schrodinger representation of the quantum mechanical harmonic oscillator. Let q and p 
be the position and momentum operators on Ti. The annihilation operator a is defined by 
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--(ojq + ip), (2.1) 



where u is the angular frequency, and its adjoint a* is the creation operator. Then the 
number operator N is defined by 

N = a ] a. (2.2) 

The number operator N has the complete orthonormal basis {|n) | n — 1, 2, . . .} of Ti. for 
which N \n) = n \n). The Hamiltonian H of the system is given by H = hu(N + |). 

In his original description of the quantized electromagnetic field, Dirac [0] postulated 
the existence of a Hermitian phase operator (f)£> such that the unitary exponential operator 
expi0£> of <j>D would appear in the polar decomposition of the annihilation operator 

a = (expi4> D )N- 1/2 . (2.3) 

The difficulty with this approach were clearly pointed out by Susskind and Glogower by 
showing that the polar decomposition of a can be realized by no unitary operators. Instead, 



they introduced the partial isometries representing the exponentials of the phase 

exp SG z0 = (N + iy^a, (2.4) 
exp 5G -i0 = a f (iV+ l)' 1/2 , (2.5) 

and the Hermitian operators representing the sine and cosine of the phase 

cosset = ^(exp SG i(p + exp SG -i(f)), (2.6) 
sin 5G = ^-.{€xp SG i<p-€xp SG -i<j)). (2.7) 

Note that the places of the carets in these Susskind-Glogower (SG) operators suggest that 
these operators are not derived by the function calculus of a certain Hermitian operator 
corresponding to 0. These operators are considered to behave well in the classical limit, 
but they fail to define well-behaved operators even for periodic functions of the phase in 
the quantum regime. Thus we cannot derive the correct statistics of the phase from these 
operators. However, it turns out that these operators give the correct mean values of the 
corresponding quantities e ^, sin<^> and cos (p. A systematic method for obtaining the correct 
statistics needs a new mathematical concept generalizing the Hermitian operators, which is 
described in the next section. 



3. Quantum estimation theory 



Denote by C(H) the algebra of bounded operators on H. In the conventional approach f36 
any observable A has a unique Hermitian operator A on TC, and any Hermitian operator A 
has a unique resolution of the identity E^(x) (x G R), called the spectral family of A. Then 
the probability distribution function of A in a (normalized) state ip G TC is given by 

Pr{A<x\?p} = (ip\E A (x)\i)). (3.1) 
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For any Borel function /, the observable f(A) has the Hermitian operator f(A) defined by 
the function calculus 

f(A)= f f(x)dE A (x), (3.2) 
and then the mean value of f(A) in state ip is 

Ex[/(A)M = Wim). (3.3) 

This framework of the statistical interpretation based on the Hermitian operator can be 
generalized to the framework based on the so-called POM, by now fairly well-known, as 
follows. 

A POM or non- orthogonal resolution of the identity on Ti is a family P(x) (x G R) of 
bounded Hermitian operators on Ti with the following conditions: 

(51) linx I ._ > _ 00 P(x) = 0, Hindoo P(x) = 1, and P(x) = lim^oo P{x -l-n" 1 ), where the limit 
is taken in the strong operator topology. 

(52) From x' < x" it follows that < P(x) < P{x') < 1. 

Note that it is usual that "POM" abbreviates the "probability operator- valued measure", 
which is more or less mathematically equivalent to the simpler notion of "non-orthogonal 
resolution of the identity" , as long as we restrict our attention to one-dimensional probability 
distributions. 

Now, we shall start with the following presupposition: corresponding to any measurable 
physical quantity X , there is a unique POM P x such that the probability distribution function 
of X in a (normalized) state ip &H is given by 

Pr[X < x\ip] = (i/j\P x (x)\ijj). (3.4) 

For an observable A in the conventional framework, the corresponding POM Pa is the 
spectral resolution of the Hermitian operator A. For any Borel function /, the operator 
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f(X) = J f dPx is defined as follows. Let dom(/(X)) be the set such that 

dom(fXx)) = {q G H\ I \f(x)\ 2 d^Px^) <oc}. 

For any rj G dom(/(X)), f{X)r] is defined by the relation 

(Z\f(X)v)= I f(x)d(Z\P x (x)\ V ), 

for all xi G Tl, where the integral is Lebesgue-Stieltjes integral. The mean value of the 
quantity f(X) in state if) G dom(/(X)) is given by 

Ex[f(XM = WiXM- (3-5) 

Then this generalizes the function calculus (|3.2|) based on a spectral resolution. Namely, 
Px is an orthogonal resolution, i.e., Px{x) is a projection for each x G R, if and only if 
X — J xdPx is a Hermitian operator on 7i such that E x = Px and f{X) = f{X) for any 
Borel function /. 

By a measurement of quantity X, we mean any experiment the outcome of which is pre- 
dicted by the probability distribution function given by ( |3.4| ). For measurability of quantity 
X with POM Px(x), it is known that for any POM Px{x) on 7i, there is another Hilbert 
space JC, a unit vector £ G /C, a unitary operator U on 7i ® JC, and a Hermitian operator A 
on /C with spectral resolution satisfying 

(il>\Px(x)\il>) = (V ® 61^(1 ® J5 A (as))Z7|V ® 0, (3-6) 

for all if) G 7^ [p7[| . This experiment consists of the following process; 1) preparation of the 
apparatus (described by JC) in state £, 2) interaction (described by U) between the object 
(described by 7i) and the apparatus, 3) measurement of the observable A (corresponding 
to Ea(x)) in the apparatus. We shall call any experiment with the above process which 
satisfies Q3.6D as a measurement of quantity X or POM Px- Thus for any quantity X with 



POM Px there is a measurement the outcome of which is predicted by ( |3.4| ). Note that this 
interpretation of the statistics of the outcome of the measurement is a natural consequence 
of the conventional postulate ( |3.1| ). Our presupposition is thus a conservative extension of 
the conventional formulation of quantum mechanics, in the sense that, if every observable 
with the Hermitian operator can be measured, so can every quantity with the POM. 

The determination of the statistics of the phase is thus reduced to the determination of 
the POM P^ corresponding to the phase <fi. This problem is solved in quantum estimation 
theory as follows. We assume that the phase <fi has values in [0, 2ir), and hence we require 
P(j,{x) = for x < and P^tt) = 1. Since the phase is canonically conjugate to the action 
in classical mechanics, the number operator is the infinitesimal generator of the phase shift 
operators e l9N in quantum mechanics. Thus the POM P^ should satisfy the relations 

e- i9 "dPJ9')e ie " = dPJ6' e 9), 

(3.7) 

< 6', 6, 6' e 6 < 2tt, 6' Q 6 = 6' - 9 (mod 2tt). 

Or equivalently, 

e-»*P*{B)e»» = PitB-*), 

(d.tt) 

B,B„g C [0,27r), B^ e = B-9 (mod 2tt), 
where P^B) = J B dP^. Any POM satisfying fl377|) is called a covariant POM 0. It is 
well known that there is no Hermitian operator such that its spectral resolution satisfies 
the above relations, but there are many solutions among general POM's. In order to select 
the optimum one, consider the following estimation problem of the phase parameter 9. Let 
us given an optical mode in a reference state if) G Tt which is supposed to interact with a 
phase shifter with unknown shift parameter 9 (0 < 9 < 2tt) so that the outgoing state is 
ipe — e t0N ip. The estimation problem is to find an experiment in state ipg which gives the best 
estimate of the parameter 9. This is equivalent to find a measurement in the state if>e the 
outcome 9 of which is the best estimate of the parameter 9. The relevance of this estimation 
problem to the determination of the POM for the phase is as follows. Suppose that the 



reference state ip were the phase eigenstate ip — \<p — 0). The outgoing state from the phase 
shifter would be the phase eigenstate ipe = \<p = 9), for which the the best estimator would 
give the estimate 9 = 9 with probability 1. Thus in this case the best estimate results from 
the measurement of the phase. Thus for a POM P to represent the phase, it is necessary that 
it is the optimum estimator of this estimation problem if the reference state approximates 
the phase eigenstate closely. Whereas we do not know what and where are phase eigenstates, 
we can reach the essentially unique solution as follows. For a POM P, the joint probability 
distribution 

d9 

p(de,d&) = {if> e \dP(§)\i/>e) 7r (3.9) 

gives naturally the joint probability distribution of the true parameter 9 and the estimate 9. 
Given an appropriate error function W(9 — 9), which gives the penalty for the case 9^9, 
the optimum estimator should minimize the average error 

/ W(9-9)p(d9,d9). (3.10) 
Jo 

Optimization problems of this type have been studied extensively in quantum estimation 
theory || [/J. The following POM P opt is the covariant POM which is the common optimum 
solution for a large class of error functions such as W(x) = 4 sin 2 | or W(x) = —5(x), where 
8(x) is the periodic 5-function: 

d9 

(n\dP opt (9)\n ! ) = e K<xn-* nl ) e i{n-n<)6 ^ ( 3 _ n ) 

where \n) (n = 0, 1, . . .) is the number basis and a n = arg(n|^)- Note that an optimum 
solution for the particular error function —5(x) is called a maximum likelihood estimator, 
and for an arbitrary error function W a Bayes estimator for W. Since this problem is not 
of the estimation of the absolute phase of the outgoing state, the optimum solutions depend 
on the phase factors a n of the reference state ip. However, this dependence only reflects 
the our optional choice of the phase eigenstate \<p = 0), and each choice of the optimum 
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POM P opt determines a unique |0 — 0) among physically equivalent alternatives. To see this 
physical equivalence, replace each number state |n) by the physically equivalent e lOLn \n) 1 and 
the same POM P opt turns to be the solution for a n — 1. Thus the particular choice of a n 
does not affect the physics. For simplicity, we choose the solution for a n — 1 (n = 0, 1, . . .) 
and determine it as the POM P^ for the phase 0, i.e., 

(n\dPJ9)\n') = e i{n - n ' )d — . (3.12) 

We call P^ as the phase POM. 

Applying (|3.4|) and ( |3.5|) to the phase POM, we obtain the statistics of the phase as 
follows. The probability distribution function of the phase in state %fj = J2 n c n\n) G Ti is 

Pr[0<^] = £c n <, J e S n-ny^ 
n,n' ® 

The mean value of the Borel function /(0) of the phase is (?/>|/(0)| , 0), where we have 

m = jT f(0)dP46) 

= £|n>/ e i{n - n ' )e m7r(n'\. 

, JO 27T 

An interesting result from this is that the SG operators coincide with the operators defined 
from the phase POM 0, p. 141], i.e., 

e5Ep SG ±z0 = e ±w dP46) = e±#, 

j o 

/■27T _ 

COS5G0 — / cos 9 dPAO) = COS 0, 
JO 

_ /-27T . „ 

sinsc0 = / sin 9 dPJ9) = sin 0. 
jo 

Thus the SG operators give the correct mean values of /(0) = e ±i<?i , sin0 and cos0, but none 
of their powers. 

The phase POM P^ gives the correct mean value for all Borel functions of the phase, but 
it gives little information about the algebraic structures of the physical quantities including 
the phase. 
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(3.13) 
(3.14) 
(3.15) 



4. The Pegg-Barnett operators 

Now we shall turn to the proposal due to Pegg and Barnett O. They start with the s- 
dimensional subspace \I/ S of Ji spanned by number states \n) with n — 0, 1, . . . , s — 1. For 
9 m = mA9 (m — 0, 1, . . . , s — 1), where A9 = 2n/s, the approximate phase state is 

\9 m )=s- 1 / 2 Y / e me -\n), (4.1) 

n=0 

and the approximate phase operator <f) s on is 



s-l 



E^IO^ml- (4-2) 
m=0 



Their intrinsic proposal is that the mean value of the quantity /(0) in state tp is the limit 

of (ip\f((i>s)\il>) as s — > oo. Then for state ^ = Sn=o c n|^) (fc < oo), we have 



A/9 

ra,n'=0 m=0 Z7F 



s-l 

D i{n— n') 



k r 2n Aft 

n,n'=0 



for any continuous function /(#) on [0, 27r]. Thus, the mean values are the same as those 
given by the phase POM P^[d9\ and we have 

lim f(l) = r f(9)P 4> (d9) = /(J), (4.3) 

where the limit is taken in the weak operator topology. In particular, the limit of their 
exponential, sine and cosine phase operators are the SG operators, i.e., 

lim exp ±i<p s = exp SG ±i<f), 
lim cos d> s = COS5G0; 
lim sin0 s = smsG(f>- 
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Therefore, the statistics of the phase obtained by Pegg and Barnett coincides with the 
statistics obtained by the phase POM, and that the limit of exponentials of the approximate 
phase operators on finite dimensional spaces is nothing but the SG exponential operators, 
as long as the limit is taken in the weak operator topology on the Hilbert space 7i. 

In the following sections, we shall develop an entirely new approach to the limit of the 
approximate phase operators. 

5. Nonstandard analysis 

It seems that certain amount of a physical quantity can be described both in quantum me- 
chanics and classical mechanics consistently. Consider a highly excited single-mode radiation 
field with the number n of quanta. Let E q the energy of this radiation field in quantum me- 
chanics with unit U q , and E c in classical mechanics with unit U c . Then, if both mechanics 
describe the same physical state, we should put E q U q = E C U C . According to quantum me- 
chanics, we have E q = % q ui q {n q + |). We can choose the unit system of quantum mechanics 
so that H q — 1, and the unit of the time is assumed to be common in both mechanics, so 
that oj q = uj c . In classical mechanics, the position-momentum uncertainty is negligibly small, 
and hence it is necessary to choose the unit system of classical mechanics such that H c is 
negligibly small. In classical mechanics, the relation E c = h c u c (n c + \) is consistent with 
the fact that the classical energy E c is a continuous variable, since the minimum increment 
Ti c oj c is negligibly small. However, in order for E c to be a finite positive number, n c should 
be infinitely large. If we do not invoke negligibly small h c nor infinitely large n c in classical 
mechanics, nothing appears to be inconsistent. However, a question arise — can we assume 
that n q is finite and yet n c is infinite? Since the number of quanta is a dimensionless quantity, 
we should put n q = n c , and hence n q should be also infinitely large. This consideration sug- 
gests that if we extend quantum mechanics to the states with infinite number of quanta, the 
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correspondence between quantum mechanics and classical mechanics becomes much more 
consistent and thorough. The rigorous mathematical description of this extension of quan- 
tum mechanics has been already possible by use of nonstandard analysis. In this section, the 
basic principles of nonstandard analysis is explained, and in the next section the macroscopic 
extension of quantum mechanics is described. 

In nonstandard analysis, we consider a large collection V of sets which contains every set 
used in the ordinary analysis. Note that V contains a set of a set, a set of a set of a set, etc. 
Any set M in V has its nonstandard extension *M. A set X which is an element of some 
*M is called an internal set. We have two basic principles in nonstandard analysis, called 
the transfer principle and the saturation principle. The transfer principle states that any 
elementary property which holds for Mi, M2, . . . , M n also holds for *Mi, *M2, . . . , *M n . A 
difficult point in nonstandard analysis is to understand the notion of an elementary property. 
The rigorous explanation of this notion needs a substantial part of mathematical logic. 
Fortunately, most of mathematical properties used in physics is elementary properties — 
a reason why the intuitive use of infinitesimal numbers has been useful in physics — and 
hence "any elementary property" can read "any property" in most of time. The saturation 
principle which we assume in this paper states that any sequence of internal sets with the 
finite intersection property has a nonempty intersection. For further detail on the basic 



framework of nonstandard analysis we shall refer to Hurd-Loeb ||21|| . Our framework is 
called an Ki-saturated bounded elementary extension of a superstructure which contains R 
in the usual terminology. 

Let *N be the nonstandard extension of the set N of natural numbers. Then *N satisfies 
all elementary properties of the natural number system by the transfer principle. An element 
of *N is called an internal natural number, and an element of *N\N is called a nonstandard 
natural number, where *N \ N is the set of elements in *N which are not in N. Every 
nonstandard natural number is greater than any (standard) natural number, and yet has 
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all elementary properties of (standard) natural numbers. An internal set enumerated by a 
nonstandard natural number with an internal one-to-one function is called a hyperfinite set. 
Every hyperfinite set is an infinite set, but has all elementary properties of finite sets. 

Let *R be the nonstandard extension of the set R of real numbers. Then, *R is a proper 
ordered field extension of R by the transfer principle. An element of *R is called a hyperreal 
number. The hyperreal number system *R is not complete nor satisfies the Archimedean 
axiom, since these properties are not "elementary properties". A hyperreal number x is 
called infinite if \x\ > n for any n G N, finite, write \x\ < oo, if there is some n G N such 
that \x\ < n, and infinitesimal if \x\ < 1/n for any n G N. For any x, y G *R, we will 
write x ~ y if \x — y\ is infinitesimal. For any finite hyperreal number x, there is a unique 
real number r such that r m x; this r is called the standard part of x and denoted by °x. 
Any function / from X to Y is extended to an internal function */ from *X to *Y. A 
sequence a n G R {n G N) is extended to an internal sequence *a u G *R {y G *N), so that 
Hindoo a n = a if and only if *a v ~ a for all v G *N \ N. A function / : dom(/) — > R with 
dom(/) C R is continuous at x G dom(/) if and only if *f(x) ~ *f(y) for all y G dom(/) 
with x ~ y. The above notions and results of infinitesimal analysis are naturally extended 
to the nonstandard extension *C of the complex number field C. 

Let (E,p) be an internal normed linear space over *C with norm p. We define the 
principal galaxy Eq of (E,p) and the principal monad Em of (E,p) as follows: 

E G = {x G E | p(x) < oo}, (5.1) 
£ M = {a; G £ | p(x) w 0}. (5.2) 

Then both and are linear spaces over C. Let E = E G /E M and p(°x) = °(p(x)) 
for x G where °i = x + Then (E,p) becomes a normed linear space over C, 

called the nonstandard hull of (E,p). By the saturation principle, it is concluded that the 
nonstandard hull {E,p) is a Banach space |2"I , p. 156]. When E has an inner product (-|-) 



14 



which determines the norm p, the nonstandard hull (E, p) is a Hilbert space with the inner 
product such that ( °x\ °y) = °{x\y) for all x, y G Eq. 

6. Macroscopic extension of the Schrodinger representation 

Let *H be the nonstandard extension of the Hilbert space H of quantum states, and *H its 
nonstandard hull. Then is a Hilbert space with inner product ( °£| °r\) = °{£\r)) for i] E 
*TCg- Let v be a nonstandard natural number and D the internal //-dimensional subspace 
of *H spanned by the hyperfmite set {\n) \ n — 0, 1, . . . , v — 1}. Then the nonstandard hull 
T> contains 7i as a closed subspace by the canonical isometric embedding V-h which maps 
^GHto °*£eT>, i.e., £ = V H £= ° % The nonstandard hulls *H and V are, thus, Hilbert 
spaces which satisfy the relations 

n c v c *n. (6.i) 

Let A be the internal algebra of internal linear operators on V. Then A is a hyperfmite 
dimensional internal *-algebra over *C. For x G A, let ||x|| be the internal uniform norm 
of x. Then, the nonstandard hull A of (A, \\ ■ ||) becomes a Banach *-algebra with norm 
|| = ||x||. Then the norm || • || satisfies the C*-condition, i.e., || ( °rc)* (°x)|| = || °x\\ 2 for 
all °x G A, and hence (A, \\ • ||) is a C*-algebra. 

Any internal operator x G Ac leaves V G and V M invariant and gives rise to a bounded 
operator n(x) on V such that ir(x) °^ = °(x£) for all £ G V G . Then the correspondence 
°x i— > 7r(x) for a; G ^ defines a faithful *-representation of the C*-algebra A on the Hilbert 
space V. 

Let Px» be the internal projection from *H onto V. Any bounded operator T onH has 
the nonstandard extension *T which is an internal bounded linear operator on *H. Denote 
the restriction of Pv*T to V by *Tx>. Then it is easy to see that *T-p G ^4g- We denote the 
operator 7i(*Tx>) on V by Tp, which is called the standard hyperfinite extension of T to P. 
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Then T-p = T on H. and ||Td|| = \\T\\. Properties of such extensions from TitoV are studied 



by Moore [5^] extensively. Now the following statement is easily established; cf. p8| , Lemma 
1.3]. 

Theorem 6.1. The mapping Sv '■ T i— > is a completely positive isometric injection 
from C-iTi) to A and the mapping E^t : it(T) i— > V^7i(T)V^ is a completely positive surjection 
from A onto C(7i). The composition £n£v is the identity map on C(TC), and the composition 
£v£h is the norm one projection from A onto £t>(£(H)). 

Let A be an observable represented by a bounded Hermitian operator A on the Hilbert 
space 7i of the Schrodinger representation. Then as described by Theorem |6.1| , the operator 
A is extended to the bounded Hermitian operator Ad on the space T>. Thus, we have a new 
quantization which associates an observable A with a Hermitian operator on T>. The above 
theorem states that this quantization leaves the mean values of observables in quantum states 
of the Schrodinger representation unchanged, i.e., 

for any bounded observable A and state ip G TC. 

Our claim is that there is a physical quantity such that its Hermitian operator on Ti does 
not exists, but its Hermitian operator on T> does exist. In the next section, we show that the 
phase is one of such physical quantities. In Section 8, we show the macroscopic character 
of states in T> orthogonal to Tl, which suggests the semiclassical nature of the observables 
which fail to have the Hermitian operator on 7i. 



16 



7. The Hermitian phase operator 

Let A9 = 2tt/u, and 9 m = mA9 for each m (m = 0, 1, . . . , v — 1). The internal phase 
eigenstate \9 m ) in V is defined by 

\9 m ) = v- 1/2 J2e m6 ™\n). (7.1) 

n=0 

Then we have 

\Qm\Qm') = $m,m'- (7-2) 

The internal phase operator <pi on T> is defined by 

v-\ 

4>I = ^2 ®m\9 m )(9 m \. (7.3) 
m=0 

Then the internal phase operator (pj has the internal spectrum {2m7T jv \m = 0,1, . . . ,v — 1} 
and hence is in *4.o Thus we have the Hermitian operator vr(0/) on X>, denoted by and 
called the Hermitian phase operator on XX Denote by A(^) the spectrum of cf> and n o (0) the 
point spectrum (eigenvalues) of </>. 

Theorem 7.1. We have A(0) = n o (0) = [0,2?r]. For eac/j 6 R (0 < < 2tt) ; 
vector °\9 m ) G X? 9 m ~ 9 is an eigenvector of <f> for the eigenvalue 9. 

For each 9 G R (0 < < 27r) and n G N, define F{6, n) to be the internal projection 

F(9,n)= ]T K)(9 m \. 

Then F(9,n) = if 9 + n" 1 < 0, and for each 9 the sequence n(F(9,n)) (n G N) is a 
monotone decreasing sequence of projections on X>. Define E^{9) to be the strong limit 



of n(F(9,n)). Then E^iff) (9 G [0,27r]) is the spectral resolution for <p fl38| , Theorem 4.1]. 
Now it is shown in the next theorem that the Hermitian phase operator has the desirable 
properties. 
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Theorem 7.2. The Hermitian phase operator <p on T> with its spectral resolution E$ 
satisfies the following conditions (P1)-(P4): 

(PI) The spectral resolution of (p is a Naimark extension of the phase POM P^, i.e., 
for any state ip G H. 

(P2) For any continuous function f , the Hermitian operator /(</>) is the limit of f((f) s ) as 
s tends to infinity in the weak operator topology of H, where 4> s is the Pegg-Barnett 
operator, i.e., 

for all states ip ETC. 
(P3) The Susskind-Glogower phase operators are given by the relations: 

exp 5G ±i0 = V-h exp±i<t>VH, 

COS5G0 = V^COS(f)V H , 

sin 5G = Vh sin(f)V n . 
(P4) For any continuous function f , 

= T,*f(Om)\(o m m\ 2 , 

m=0 

for all ip G T>g- 

As consequences from (PI), we obtain the statistics of the phase only from the ordinary 
quantum rules with the Hermitian operator <p. Let ip G Ti be a normalized state of the system. 
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Then ip is in the domain of by relation ( |6.1|) . The probability distribution function of the 
phase in the state ip is given by 

Pr[0<^] = (^(W>- 
For any Borel function /, the physical quantity /(</>) has the Hermitian operator 

/(#)= f(9)dE <j> (9), 
Jo 

by the function calculus, and the mean value of f(<p) in the state ip is given by 

ex[/(0)^] = w$m. 

8. Macroscopic states 

We have extended the Schrodinger representation on 7i to the hyperfmite dimensional space 
T>. The following theorem shows that the states in T> TC can be interpreted naturally as 
the classical limits of the quantum mechanical states in 7i. 

Theorem 8.1. Let T G C(7i). Suppose that (n\T\n') (n, n' G N) is a Cauchy sequence 
in n andn 1 . Then for any nonstandard k, k' e *N with \k), \k') 6 D the standard hyperfinite 
extension T v E A of T satisfies the relation 

{k\T v \k') = lim (n\T\ri). (8.1) 



Remark. Suppose that the nonstandard universe is constructed by a bounded ultra- 
power of a superstructure based on R, with the index set I = N and a free ultrafilter IA. 
Then any nonstandard number k G *N \ N is represented by a sequence s(i) (i G N) of 
natural numbers in such a way that two sequences s(i) and s'(i) represents the same non- 
standard number k if and only if {i G N| s(i) = s'(i)} G U. Let m be a nonstandard number 
corresponding to a sequence s(i). Let a(n) be a bounded sequence of complex numbers. 
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Then the standard part of *a(m) coincides with the ultralimit of the subsequence a(s(i)) of 

a (*a(m)) = lim a(s(i)). (8.2) 

i— >W 

Thus, for any T G C(TC) and fc, fc' G *N \ N, we have, 

(k\T v \k')= lim(s(i)\T\s'(j)), (8.3) 

provided k, kl are represented by sequences s(i), s'(j). Thus the matrix element (k\Tx>\k') is 
the ultralimit of a subsequence of (n\T\n') (n, n' G N), even if it is not a Cauchy sequence. 

Let A be an observable represented by a bounded Hermitian operator A on the Hilbert 
space TC of the Schrodinger representation. Since Ad is an extension, the matrix element 
(n\Ax>\n') for standard n, n' G N is the same as that of the operator A, i.e., {n\A-o\n') = 



(n\A\n'). From Theorem |S.1| , the matrix element (n\A-p\n') for nonstandard n, n' has also 
a clear physical interpretation, i.e., it is the classical limit of the matrix elements (n|A|n'). 
Thus, the states in TC \ T> are naturally considered as the classical limits of quantum states. 
We shall call TC the purely microscopic part of T>, and T> TC the purely macroscopic part 
of V. The space T> contains microscopic states in TC, macroscopic states in TC f>, and 
superpositions of those states in T>. Thus, it is appropriate to call the extension TC —*T> with 
the extensions A — > Ap for observables A as the macroscopic extension of the Schrodinger 
representation. The greatest significance is, of course, that we can associate the phase with 
a Hermitian operator on the macroscopic extension. 

9. Conclusions 

The claim, "We have recently shown that a Hermitian optical phase operator exists. This 
result contradicts the well established belief that no such operator can be constructed." due 
to Barnett and Pegg || is not appropriate. Instead, they construct a sequence of Hermitian 
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operators (f) s on finite s-dimensional spaces such that the statistics obtained by the operators 
4> s approaches to the statistics of the phase as s tends to infinity. Their statistics of the phase 
obtained by this limit process is consistent with the one obtained from the phase POM. 

In this paper, we give a demonstration of the claim that the correct statistics of the phase 
is the one given by the phase POM as well as the Pegg-Barnett limit process, and construct a 
single Hermitian operator <p on an infinite dimensional Hilbert space which gives the statistics 
of the phase by the ordinary quantum rules. Furthermore, the Hermitian operator is 
infinitesimally close to the internal phase operator (pj, that is, the statistical prediction given 
by (f> is obtained from <fii with only infinitesimal difference. The remarkable character of 
the internal phase operator <fii is that this operator shares all elementary properties of the 
Pegg-Barnett operators on finite dimensional spaces. This follows from the transfer principle 
in nonstandard analysis. 

Thus the statistics of the phase variable is obtained by the following procedure. Fix a 
nonstandard natural number v and consider the //-dimensional internal linear space T>. The 
internal phase operator <pi is well-defined on T>. Recall that all number states and all phase 
states are contained in T> (cf. Theorem [7.1| ). Every computation involving 0/ is formally the 
same as the computation involving the Pegg-Barnet operator, but gives the correct statistics 
with only infinitesimal error. The last step is, if desire, to remove infinitesimals from the 
result of the computation. 



Theorem [7.2| states that this procedure is equivalent to the computation with the usual 
quantum rules involving the Hermitian phase operator on the Hilbert space T>, which 
extends the Hilbert space 7i of quantum states by adding macroscopic states. This result 
does not contradict the well established belief that no such operator can be constructed on 
the Hilbert space of quantum states. Instead, we have shown that the phase operator can be 
constructed on the Hilbert space of the quantum states plus the macroscopic states. 
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A. Appendices 

A.l. Proof of Theorem 7.1 

Let 9 G [0, 2n] and 9 m ps 0. Then obviously, 

and hence °\9 m ) is an eigenvector of <fi corresponding to eigenvalue 9. Thus we have [0, 2n] C 
n o (0). Since < 9 m < 2n, we have < {ip\(f>i\ip) < 2n for any unit vector ip E T>, and hence 

It follows that < < 2?rl, so that A(0) C [0,2tt]. Therefore, A(0) = [0, 2tt] = n o (0). This 
completes the proof. 

A.2. Proof of Theorem 7.2 

Let n, ml G N. We have 

(n|^(0)|ra') = lim (n\F(9,k)\n') 

k— >oo 

= lim °(n\F(9, k)\n) 

k—*oo 

= lim V e^- n '^ — 







6 e i(n-n>)5 

2n' 
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From ( |3.12j ), this concludes condition (PI). Let a s (s G N) be a sequence such that a s = 
(n\f(<j) s )\n'), and *a s (s G *N) its nonstandard extension. Then we have lim^oo a s ~ *a u , 
and hence 

Hm(n|/(0 s )|n') « (n| */(0/)|n') 

f — 1 A/3 

Jo 2n 
= (n\f$)\ri). 

Since the first and the last term is standard, they must be the same, i.e., 

\im(n\f(4 >s )\n') = (n\f(4 > )\n'). 

Then (P2) follows easily, To prove (P3), by the obvious relations, it suffices to show the first 
equation. We have 

(n|e-*V) = {n\dE^ff)\n') 

f 2n i(n-»'-l)0 ^_ 

Jo 2tt 

— On,n'+l: 

whence the desired equation is obtained from the relation 

(n\(N + l)-*a\n') = (n\n' + l) 

Thus (P3) follows easily. (P4) is obvious, and the proof is completed. 
A. 3. Proof of Theorem 8.1 

Let a n>n i = (n\T\n') for n, n' G N, and L = rim Tl)Jl /_K 30 a n>n i. Let *a mjm > ( m, m' G *N) 
be the nonstandard extension of the sequence a n<n i (n, n' G N). Then, L « a mjm / for all 
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m, ml E *N\N. Let k, k' e {0, 1, . . . , v- 1}\ *N. By transfer principle, *a fc , fc ' = *T|fc'), 
and hence L « *T|A;'). It follows that L = °(k\*T\k') = (k\T D \k'}. This completes the 
proof. 
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